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The binding energies of the ground states and several excited states related to single-particle and 
-hole states in nuclei around ^^O are calculated taking charge dependence into account. Effective 
interactions on the particle basis are constructed from modern charge-dependent nucleon-nucleon 
interactions and the Coulomb force within the framework of the unitary-model-operator approach. 
Single-particle (-hole) energies are obtained from the energy differences of the binding energies 
between a particle (hole) state in ^^O or ^^F (^^N or ^^O) and the ground state of ^^O. The resultant 
spin-orbit splittings are small for the hole state and large for the particle state in comparison with 
the experimental values though the differences between the experimental and calculated values are 
not very large. The charge dependence of the calculated single-particle energies for the ground states 
are in good agreement with the experimental values. Furthermore, the Thomas-Ehrman shift due 
to the Coulomb force for the lsi/2 states in ^^O and ^^F can be observed. 

PACS numbers: 21.30.Fe, 21.60.-n, 21.10.Pc, 27.20.-|-n 
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I. INTRODUCTION 

The single-particle level is one of the fundamental 
structures in nuclei. Important physical quantities such 
as the spin-orbit splittings and the magic numbers are 
characterized by the single-particle level. Recently, it has 
been argued that some magic numbers disappear and new 
magic numbers arise in nuclei near the drip lines P, . 
When we calculate the energies of single-particle levels in 
neutron- or proton-rich nuclei, it would be desirable that 
the calculation formalism is based on the particle basis. 
Advantages of the particle-basis formalism are that the 
Coulomb force can be treated accurately for the proton- 
proton channel and effects of charge dependence in re- 
alistic nuclear forces are taken into account in structure 
calculations. In the particle-basis formalism, one can ob- 
tain the energy differences between proton and neutron 
levels for not only N = Z nuclei but also neutron- or 
proton-rich nuclei in the same manner. 

The calculation of single-particle energies starting with 
a nucleon-nucleon force in free space is a fundamental 
problem in theoretical nuclear physics. There have been 
many attempts to understand the structure of single- 
particle levels as well as other ground-state properties 
in nuclei 0,0, IE IE 0- I'^ ^^^h calculations, we need a 
many-body theory that leads to an effective interaction in 
a restricted model space for a nucleus in many cases. For 
this purpose, the G matrix has been widely used as a ba- 
sic ingredient in performing structure calculations |E ■ 
A recent study of the doubly closed-shell nucleus ^^O 
using the G matrices constructed from modern nucleon- 
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nucleon interactions can be seen in Ref. [I^- Lately, 
as an alternative of the G matrix, a low-momentum 
potential Mow-fc has been constructed from a realis- 
tic nucleon-nucleon interaction using a renormalization 
group technique or conventional effective interaction the- 
ory by Bogner et al The application of Viow-fe to the 
calculation of ground-state properties of the closed-shell 
nuclei ^^O and ^"Ca has been done in Ref. [T^ . 

As one of the methods for solving nuclear many-body 
problems, we have developed the unitary-model-operator 
approach (UMOA) 1^. An energy-independent and 
Hermitian effective interaction is derived through a uni- 
tary transformation of an original Hamiltonian. Nuclear 
ground-state properties, such as the ground-state energy, 
charge radius, and single-particle energy have been cal- 
culated for ^''O 14J and ^''Ca ^1^]. We have learned that 
spin-orbit splittings for hole states are enlarged by tak- 
ing second-order diagrams into account. Furthermore, 
the UMOA has been developed for the structure calcula- 
tion of A hypernuclei and applied to j^O, ]^0, a nd 1 ^ Ca 
using hyperon- nucleon interactions in free space |l6t Il7l | . 
Differences in the properties of modern hyperon-nucleon 
interactions have been disclosed in the structure calcula- 
tion. 

Recently, we have extended the formulation of the 
UMOA from the isospin basis to the particle one for the 
purpose of the charge-dependent calculation. To con- 
firm the validity of the calculation method based on the 
particle basis, in this paper, we apply this method to 
^^O and its neighboring nuclei ^^N, ^^O, ^'''O, and ^'''F. 
Binding energies of these nuclei are calculated for the 
ground states, and excited states which have the single- 
particle or single- hole structure as the main component. 
The single-particle energy in the neighboring nuclei is 
given as the relative energy between a single-particle (- 



2 



hole) state in the neighboring nuclei and the ground state 
of ^^O. As for the single-particle (-hole) state, the exci- 
tation up to two-particle one-hole (one-particle two- hole) 
from the unperturbed ground state of ^^O are taken into 
account. 

Four high-precision nucleon-nucleon interactions rep- 
resented in momentum space are employed, namely, the 
Nijmegen 93 (Nijm 93) [3, Nijm I the charge- 

dependent Bonn (CD Bonn) and the next-to-next- 
to-next-to-leading order (N'^LO) potential "21] based on 
chiral perturbation theory |2^ |2^ which has recently 
been constructed by Entem and Machleidt. In these 
potentials, effects of charge dependence are taken into 
account. The first three potentials are based on meson- 
exchange models in which several kinds of meson are in- 
corporated. On the other hand, the essential degrees of 
freedom of the mesons in the N'^LO potential are only 
for the pions. Therefore, the N'^LO potential is con- 
structed in a low-momentum region compared to the 
meson-exchange potentials which have heavier mesons. 
However, the N'^LO potential has the high accuracy to 
reproduce the nucleon-nucleon data below i?iab = 290 
MeV, and thus the N'^LO potential as well as other high- 
precision nucleon-nucleon interactions can be used in nu- 
clear structure calculations. 

This paper is organized as follows. In Sec. ^ the 
methods for deriving effective interactions and perform- 
ing structure calculations are given. In Sec. IIIIL calcu- 
lated results for ^^O and its neighboring nuclei using the 
four realistic interactions are presented. Finally, we sum- 
marize the present work in Sec. IIVI 



II. METHOD OF CALCULATION 

In the UMOA, the Hamiltonian to be considered is 
given by a cluster expansion of a unitarily transformed 
Hamiltonian. In the previous works |l3lll4|. many-body 
correlations up to three-body cluster terms have been 
evaluated. It has been confirmed that the cluster expan- 
sion in the numerical calculation for ^^O shows the good 
convergence at the three-body cluster level. We may say 
that since we consider the iV ~ Z nuclei around ^^0 in 
the present study, the three-body cluster terms do not 
have a significant contribution to the energy difference 
between the single-particle (-hole) levels of the proton 
and neutron. Therefore, in the present calculation, we 
neglect the three-body cluster terms for simplicity. The 
evaluation of the three-body cluster terms based on the 
particle basis is a further challenge which should be ac- 
complished for a deeper understanding of nuclei. 

In the following subsections, we present a general 
framework for deriving an effective interaction and a 
practical method for the structure calculation in the 
present study. 



A. Derivation of effective interaction in the P and 
Q spaces 

In the usual sense of effective interaction theory, an ef- 
fective interaction is defined in a low-momentum model 
space {P space). However, in general, one can also derive 
an effective interaction in the complement (Q space) of 
the P space by making the decoupling condition for the 
effective interaction v as QvP ~ 0. Note that the pro- 
jection operators P and Q satisfy the usual relations as 
p + Q ^ 1, P^ = P, ^ Q, and PQ ^ QP = 0. We 
here present a general framework for deriving an two- 
body effective interaction of Hermitian type for a two- 
body system. 

The two-body effective interaction vi2 of Hermitian 
type is written as 

ii2 = U~'^{hQ + vi2)U -ho, (1) 

where vi2 is the bare two-body interaction and ho is the 
one-body part of the two-body system which consists of 
the kinetic energy ti (12) and, if necessary, the single- 
particle potential ui {U2) as Hq — ti + ui + t2 + U2. The 
operator U for the unitary transformation of Hq + V12 can 
be written as [2^ 

[/ = (1 +CJ - wt)(i +t^^t +^t^)-l/2 (2) 

by introducing the operator uj satisfying lo ~ QluP and 
thus oj^ = oj^ — 0. The above expression of U agrees 
with the block form using the projection operators P and 
Q of Okubo 25] given by 




P(l+t^t^)-i/2p -Ftjt(l+tjcjt)-i/2g \ 

Qw(l+tjt^)-l/2p Q(l+tJt^t)-l/2Q J- 



(3) 

We should note here that the operator U is also expressed 
as 

U = (4) 

where S is anti-Hermitian and given under the restrictive 
conditions PSP = QSQ = by 

S — arctanh(cj — uj^). (5) 

In order to obtain the matrix elements of cj, we first 
solve exactly the two-body eigenvalue equation as 

(/io + wi2)|$fc) = (6) 

With the eigenvector |$fe), the matrix elements of oj on 
the basis states \p) in the P space and \q) in the Q space 
can be determined as 

d 

(g|c.b) = 5](g|Q|ci>fe)(0,b), (7) 

fe=i 

where d is the dimension of the P space, and {(j>k\ is 
the biorthogonal state of \(j>k) = P\^k), which means 
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the matrix inversion = [(p'|0fc)] ^ and satisfies 

J2p{'Pk\p){p\M = Sk,k' and J2k(p'\M{Mp) = ^P,P'- It 
should be noted that the set of eigenstates {\^k),k = 
1,2,- ■ -,(!} is selected so that they have the largest P- 
space overlaps among all the eigenstates in Eq. 

Then, in order to obtain the matrix elements of U, we 
introduce the eigenvalue equation for w^cj in the P space 
as 



uj^uj\ak) = /ifc|a/c)- 



(8) 



Using the solutions to the above equation, we define the 
ket vector as 



Wk) = —oj\ak), 



which is also written as 



{qWk) = —^{q\u}\p){p\ak). 



(9) 



(10) 



Using Eqs. (|5|)- (|10|l . we obtain the matrix elements of the 
unitary-transformation operator U in Eq. ^ as 

{p'\U\p) = (p'|(l+C.tc^)-l/2|p) 
d 

= J2il + nl)-'/'{p'\ak){ak\p), (11) 



{q\U\p) = {q\Lu{l+u;^Lur'^'\p) 
d 

= J2i^ + ^,l)-'/'^,k{qWk){ak\p), (12) 



k=l 



{p\U\q) = -{p\LU^l + LULO^)-'^'\q) 
d 

= -J2i^ + ^ilr'^'^ik{p\ak){l^k\q), (13) 



k=l 



and 



{q'\U\q) = {q'\{l+u;J)-'^'\q) 



fe=l 

(14) 

Thus, the matrix elements of the effective interaction -512 
in Eq. can be written as 

k,l 

(15) 

where \k), and \l) denote the basis states in the 

P + Q space. 

The above formulation is employed for deriving the ef- 
fective interaction in the present study. Here we note that 
since we treat a many-body system, the single-particle 
potential ui (112) in ho for both particle and hole states 
is introduced to obtain a good unperturbed energy. In 
the following, a procedure for determining the effective 
interaction and the single-particle potential is given. 



B. Two-step method for the calculation of effective 
interaction 



In nuclear many-body problems, how to determine 
the single-particle potential for particle (unoccupied) 
states as well as hole (occupied) states is important in 
connection with the evaluation of many-body correla- 
tions [2^, 113, m, I23. In our calculations, the single- 
particle potential, which is determined self-consistently 
with the two-body effective interaction, is calculated up 
to a sufficiently high-momentum region. In general, this 
choice of the single-particle potential leads to a deeper 
binding of the ground-state energy of a nucleus in the 
lowest order. Then, effects of the many-body correla- 
tions of higher order become smaller than the choice of 
only the kinetic energy for the particle state. This trend 
would be favorable when the evaluation of many-body 
correction terms has to be limited in the actual calcula- 
tion. 

In our earlier calculations, the effective interaction was 
derived by a three-step procedure with some approxima- 
tions to take account of single-particle potentials up to 
a high-momentum region. In the present work, however, 
we adopt a two-step procedure and approximation meth- 
ods are refined, because the performance of the computer 
has been greatly improved and some approximations in 
the previous works are not needed at present. In the 
following, we shall give the two-step procedure for the 
numerical calculation. 



1. First-step calculation 

In this work, we employ the harmonic-oscillator (h.o.) 
wave functions as the basis states. Two-nucleon states 
for Z = nn, np, pp channels consisting of the product of 
the h.o. states are given by 



\aP)z = \nJajama,nblbjbmb)i 



(16) 



The model space P^^'' and its complement Q^'' composed 
of the two-nucleon states for the Z channel are defined 
with a boundary number pi as 



[ Q^'' otherwise. 



(17) 



which is also illustrated only for the np channel in Fig. ^ 
The nn and pp channels are considered similarly in the 
actual calculation. The value of pi is taken as large as 
possible so that the calculated results do not depend 
on this value. The pi dependence of calculated results 
will be investigated in Sec. IIIII The symbols p„ and 
Pp in Fig. n stand for the uppermost occupied states of 
the neutron and proton, respectively, and in the present 
case of ^^O, p„ and pp are the Opi/2 orbits. The Qxi 
and Qx2 spaces defined with pi, /?„, pp, and px in the 
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FIG. 1: Model space Pnp and its compliment Qnp for the 
np channel in the first-step calculation. 



Qnp space should be excluded due to the Pauli princi- 
ple when we calculate matrix elements of the bare two- 
body interaction. The value of px is determined so that 
the Pauli principle from the states in the Qxi and Qx2 
spaces can well be taken into account, and taken as 
px = 2na + Za + 2nb + Zfc = 20 in the present study. 

It is noted that, in this first step, the effective inter- 
action is constructed using the relative and cm. states 
of the h.o. wave functions. Since we consider a huge 
Hilbert space, it is very difficult to use the basis states 
composed of the product of the single-particle h.o. states 
in the model and complementary spaces. In the follow- 
ing, we shall give a practical method for calculating the 
effective interaction and the single-particle potential. In 
order to derive the two-body effective interaction for each 
oi Z = nn, np, pp channels we rewrite Eq. ^ in terms 
of the relative and cm. states as 



Ek\k-l{l')SJr,NL), 



)]\k-l{l')SJr,NL)z 
(18) 



where I {I') and S are the orbital angular momentum 
and spin of a relative state, and is the total angular 
momentum given by Jr = 1 -f S. The letter k means 
an additional quantum number specifying an eigenstate. 
The terms and tc.m. are the kinetic energies of the 
relative and cm. motions, respectively, and Vi2 is the 
bare interaction. The operator projects two-body 

states on to the Q^^^ space, but Pauli-forbidden two-body 
states in the Qxi and Qx2 spaces are excluded. The sum 
of two single-particle potentials in the relative and cm. 
states is denoted by ui2{N, L). We assume, in the present 
study, that the matrix elements in solving Eq. H18|) are 



diagonal in each of the cm. quantum numbers TV and 
L. Thus, the resultant effective interaction becomes also 
diagonal in the cm. quantum numbers. 

The matrix elements of ui2(iV, L) can be written under 
the angle-average approximation [l^ | ' 



as 



{nlSJr, NL\ui2{N, L)\n'l'SJr, NL) z 

(„1)A+A' 



E 
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jb } { lb 
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jb } 
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J j[x' 
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J J 



[L] 



W{LIJS] XJr)W{Ll'JS; X'Jr) 



X {nlN LX\nalanblbX) {n I' N LX' \n'^lani,li,X') 

X{{nala3a\ui\^\n'jaja) + {UaLj alug^ KLja)) , (19) 

where [x] = 2a; + 1 and J is the total angular momentum 
for two single-particle h.o. states given by J = ja + jb- 
The coefficients {• • •}, W{- ■ •), and {nl - ■■ \nala ■ •■) denote 
the Wigner 9-j symbols, the Racah coefficients, and the 
h.o. transformation brackets, respectively. Note that, 
as for the nn and pp channels, the calculation should be 
done only for I + S — even. The quantities ui\^ and ui^J 
represent the single-particle potentials of the neutron Un"^ 
or proton Mp^^ in the first-step calculation, depending on 
Z = nn, np, pp channels. The single-particle potentials 
Un'^ and Up^^ are calculated self-consistently with the two- 
body effective interaction, which will be shown later. 

The operator can be written under the angle- 

average approximation as 



E 

nlNLSJr 
pi<2n+l+2N+L 



z{n,l,N,L,S,Jr) 
x\nlSJr,NL){nlSJr,NL\, (20) 



where 



9z{n,l,N,L,S,Jr) 

E 



(-1)^+^7: 



abXX' J 
pi<2na-\-la-\-2nh-\-lb<P^ 



( la 


1 

! 




la \ 


ja 1 




2 


III 


k 2 


jb } 


I A 


S 




. A' s 


J J 



[x][x'][jM[sM 

[L] 



W{LIJS; XJr)W{LlJS; X'Jr) 



X {nlN LX\nala'nhlbX) {nlN LX' \nala'nblbX') 



with 



Iz = 



2 for Z = nn or pp, 
1 for Z = np. 



(21) 



(22) 
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Note that, as for the nn and pp channels, the calculation 
should be done only for Z + S* = even. The letter a (b) for 
the summation in Eq. (|21|l means a set of the quantum 
numbers a = {na,la, ja, z = n or p} of a single-particle 
h.o. state. The conditions of the summation of single- 
particle states a and b for the nn and pp channels are 
{a < pm b > pn} and {a < Pp, b > pp}, respectively. As 
for the np channel, {a < /?„, b > pp} or {a > pn, b < pp}. 
Here for example, the notation {a < p„, 6 > pp} for 
the np channel means that the summation is done for 
occupied states of the neutron and unoccupied states of 
the proton. 

It should be noted that Eq. H18|) is solved exactly 
by diagonalizing the matrix elements of several hundred 
coordinate-space h.o. basis states for each channel on 
the assumption of the diagonal cm. quantum num- 
bers. If we employ a bare interaction in momentum- 
space representation, the Fourier transformation for the 
h.o. wave function is needed in calculating the ma- 
trix elements of the bare interaction. Using the eigen- 
vector \k;l{l')SJr,NL)z, the operator ui in Eq. ^ 
can be written in terms of relative and cm. states. 
Then, the matrices of the effective interaction v[^2 i^i 
Eq. (|15|) are obtained in the relative and cm. states as 
{72lSJr\vi2iN,L)\n'l'SJr)z through Eqs. Note 
that we do not need the Q-space effective interaction in 
the first-step calculation if we take a sufficiently large 
model space. 

The transformation of the effective interaction in the 
relative and cm. states into the one in the shell-model 
states can be performed straightforwardly as 

{ab\v[^^\cd)j,z 

NLW 

xV[jam[jc][jd]{m'][S][Jr] 

{la 2 Ja I I 2 ■^'^ 1 

lb I jb } { Id I jd } 
\ S J } [ X' S J } 

xW{LlJS; XJr)W{Ll'JS; X' Jr) 

X {nlN LX\nalanblbX) {n I' N LX' \nclcndldX') 

x{nlSJr\v\^2\N,L)\n'l'SJr)z, (23) 

where 

, . J 1 + {-iy+^ for Z = nn or pp, , 

which is required for the antisymmetrization of the ma- 
trix elements in the shell-model states. Note that l + S — 
even for the nn and pp channels. 

The single-particle potentials u^z^ and u^z^ in Eq. lfTO|l 
are determined self-consistently with the two-body effec- 



tive interaction Wx2 j which is written as 

m:occupicd 

x{am\v''^^\a'm) j^z (25) 
for the neutron, and 

J,Z—pp.np ^ 'V ' 

m:occupicd 

x{ma\v[2^\ma')j^z (26) 

for the proton. 

The procedure for the self-consistent calculation is as 
follows. First, we input initial values of ui\^ and ui^'^ in 
Eq. H19|l . and solve the eigenvalue equation in Eq. 118|) for 
each oi Z — nn, np, pp channels. Through Eqs. (|7I)- (|15|I . 
the effective interaction in the form of the reduced ma- 
trix element is determined. Then, the new single-particle 
potentials are calculated through Eqs. H23|l - H26|) . These 
new values of the single-particle potentials are used in 
Eq. (|19|l . and the iterative calculation is performed until 
the calculated results converge. 

We remark here that one of the practical methods of 
the structure calculations using the present effective in- 
teraction would be the shell-model diagonalization. How- 
ever, the application of such a calculation may be limited 
only to light nuclei, because we must take account of 
many single-particle states in the model space and the 
dimension of the matrices to be diagonalized becomes 
very huge. Since we intend to obtain only the energies 
of the ground state of the closed-shell nucleus and the 
single-particle (-hole) states in its neighboring nuclei, we 
proceed to the next step for a more practical calcula- 
tion. In the second-step calculation, the effective inter- 
action determined in the first-step calculation is unitarily 
transformed again so as to vanish the matrix elements for 
two-particle two-hole {2p2h) excitation. This is an essen- 
tial point of the UMOA. By virtue of this, a number of 
many-body correlations with the vertices of the effective 
interaction are reduced compared to the usual linked- 
cluster expansion with the G matrix. In the UMOA, such 
many-body correlations can be evaluated in a cluster ex- 
pansion of the unitarily transformed Hamiltonian with 
the vertices of S in Eq. (0), the one-body Hamiltonian, 
and the two-body effective interaction. 

2. Second-step calculation 

Using the two-body effective interaction -D-j^ deter- 
mined in the first-step calculation, we consider the in- 
ternal Hamiltonian as 

Hint =J2^^ + Y. -Tc.m., (27) 

i i<j 



6 



where Tc.m. is the kinetic energy of the cm. motion. 
In this second step, the calculations are performed us- 
ing the basis states of the product of the single-particle 
h.o. states. In order to remove spurious cm. states, we 
add the cm. Hamiltonian Hem. so as to constrain the 
ground-state cm. motion in the h.o. potential with the 
frequency as 



2 



(28) 



The h.o. potential Uc.m. can be written with the mass 
number A and the nucleon mass m as 



A-2 



A{A-1) 



, (29) 



where Xij — 



i(2m)r!2R2^. and x^j = \{^)n''vl. The 
definitions of the coordinates are R = ^X^i'^ii ^ 
1 (ri -I- Tj ) , and r y = — r j . We assume that the nucleon 
mass is the mean value of the neutron and proton. As 
for the value of /3c. m. in Eq- l|28|l . in the present study, 
we simply take as 0r..rn = 1 which could be acceptable as 
discussed in Refs. |3l|,|32|- Thus, the Hamiltonian to be 
considered in the second-step calculation becomes 



H = H. 



where 



A 



A{A - 1) 



(30) 



(31) 



Note that the above two-body interaction is A dependent. 

The central aim of the present study is to calculate the 
binding energies of the ground-state of ^^O and its neigh- 
boring nuclei, and to obtain single-particle energies using 
the Hamiltonian in Eq. (|30|l . To accomplish this with- 
out performing the full shell-model diagonalization, we 
proceed to the decoupling calculation again. The model 
space in the first-step calculation in Fig. ^ is separated 
as shown in Fig. [51 The model space and its compliment 

(2) (2) 

for the np channel are denoted by Pnp and Qnp , respec- 
tively. It should be noted that we solve the P-space and 
Q-space problems on an equal footing in the second-step 
calculation, using the effective interaction determined in 
the first-step calculation which has already incorporated 
the effect of the short-range correlation of the bare inter- 
action. The Pxi and Px2 spaces are the Pauli-blocked 
spaces in the second-step calculation. 

It would be worthy to mention the property of the ef- 
fective interaction to be determined in the second-step 
calculation. By taking the model and complementary 



P 




i(2) 



n 



FIG. 2: Models space Pnp and its compliment Q^p' for the 
channel in the second-step calculation. 



np 



spaces as shown in Fig. 13 the resultant effective interac- 
tion which is determined from the decoupling con- 
dition Q'^z^ Pz''^ — for Z — nn, np, pp has no ver- 
tices which induce 2p2h excitation. This is analogous 
to the Hartree-Fock (HF) condition which means that 
an original Hamiltonian is transformed so as to vanish 
the vertices of Iplh excitation. Although the vertices 
of the one-body non-diagonal matrix elements remain in 
determining the effective interaction, these non-diagonal 
matrix elements are diagonalized at the end of the calcu- 
lation. 

The eigenvalue equation for the Z channel in the 
second-step calculation which corresponds to Eq. lO can 
be written as 



{t 



,(2)- 



,(2)- 



(32) 

where \ '^k)j^,z represents a two-body eigenstate in terms 
of the basis states of the product of the single-particle 
h.o. states with a good total angular momentum and 
parity for the Z channel. We solve the above eigenvalue 
equation exactly by diagonalizing the matrix elements in 

(2) (2) 

the full space P^ + , and then obtain the matrix 
elements of U in this full space through Eqs. (|7l)- ((Tl|l . 
In addition, the matrix elements of U for the Pxi and 
Px2 spaces are given by 



{x'\U\x) = 5x,x' 



(33) 



and 



{p\U\x) = {q\U\x) = {x\U\p) = {x\U\q) - 0, (34) 
where \x), \p), and \q) are the basis states in the Pxi and 

(2) (2) 

Px2, Pz I ^-nd Qz spaces, respectively. 

The calculation procedure in the second step is as fol- 
lows. We first solve exactly Eq. H32|l by the diagonaliza- 



12\ 

tion. As the initial values of u\' and u 



,(2) 



Eq. 1221, 
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2 4 6 

number of iteration 

FIG. 3: Convergence of the unperturbed ground-state energy 
for the iterative calculation for ^^O in the second-step calcu- 
lation for pi = 12 and HQ = 14 MeV. The CD-Bonn potential 
is employed. 



we use the single-particle potentials determined in the 
first-step calculation. Through Eqs. lf7|l- 115|l . the efFec- 
tive interaction wj^j this second step in all the Pxi, 



(2) 



(2) 

and spaces is determined. Then, the 

■p 



X2, ^ z 

single-particle potentials Un' and u)^' are calculated in 
Eqs. using the efiFective interaction instead 

of {5^2^ determined in the first step, and the self-consistent 
calculation is performed iteratively until the calculated 
results converge. 

As a typical example of the convergence of the self- 
consistent calculation, in Fig. |21 we show the results 
of the unperturbed ground-state energy of ^^O in the 
second-step calculation with increasing number of iter- 
ation for pi = 12 with the CD-Bonn potential. The 
Coulomb interaction is included in the calculation. The 



unperturbed ground-state energy Eg,s. 
closed-shell nucleus is given by 



(unp.) 



of the doubly 



pn(unp.) 
g-s. 



(2jm + 1) ({m\t,\m) + ^(m|42)|^) j 

—n. n ^ 



z—n,_ 
m:occupicd 



(35) 



where \m) denotes a h.o. single-particle state for the 
hole state with a total angular momentum jV™- We also 
express the unperturbed single-particle energy iji.™'' -' for 
a state \a) as 

= {a\t,\a) + (alui^) \a) for z ^ n, p. (36) 



Note that E^g'!^^-' and E^^'p.''-' are implicitly A dependent 
due to the property of V-^p{A) in Eq. ||1JT1) . We see that 



^(unp.; 



the results in Fig. O converge when the number of inter- 
action is larger than 4. 



C. Diagonalization of the transformed Hamiltonian 

The transformed Hamiltonian determined in the 
second-step calculation does not contain the interaction 
which induces 2p2h excitation. However, there remain 
some terms inducing Iplh excitation in the one-body 
Hamiltonian, and coupling terms in the two-body inter- 
action between Ih and lp2h states for occupied states, 
and between Ip and 2plh states for unoccupied states. 
The transformed Hamiltonian to be diagonalized consists 
of the kinetic and single-particle potential parts, and the 
two-body effective interaction determined in the second- 
step calculation. As for the closed-shell nucleus, we di- 
agonalize the transformed Hamiltonian with the shell- 
model basis states, taking into account Iplh excitation 
from the unperturbed ground state. We denote the en- 
ergy shift from the unperturbed energy obtained by the 
diagonahzation by Eipih. As for the closed-shell nucleus 
plus one-particle (one-hole) system, the shell-model basis 
states are composed of the Ip and 2p\h states (l/i and 
lp2h states). The energy shift from the unperturbed 
energy obtained by the diagonalization is expressed by 
E^pih {Eip2h)- The binding energies BE for these sys- 
tems are given as follows. 



-BEC'^O) 



7^(unp.) 
g-s. 



E 



Iplh, 



and 



-i?i?("0,"F) = -Ki?2pih, 



-i?i?(i50,i5N) = 4:^,"P)-Hi?ip2h. 



(37) 



(38) 



(39) 



Thus, the single-particle energies Es.p. for the particle 
and hole states are written, respectively, as 



£;..p. (^^O, ^^F) = BErO) - BE{"0, "F) (40) 



and 



Es.p.rO, ^^N) = BEC^'O, '"N) - BErO). (41) 

In the following section, we shall present the calculated 
results of the energies using Eqs. (|35|I - H41() with some 
discussions. 



III. RESULTS AND DISCUSSION 

In the present study, the number of the h.o. wave func- 
tions which are used as the basis states is finite, and some 
approximations are made. Therefore, the calculated re- 
sults have the dependences on the h.o. energy hfl and 
the value of pi which specifies the model space in the 
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first-step calculation. In the following subsections, some 
calculated results are shown with the hfl and pi depen- 
dences. However, we search for optimal values of HQ and 
values of pi for which the calculated results almost con- 
verge to obtain the final results. 

In order to clarify differences in the properties of 
modern nuclcon-nuclcon interactions, four interactions 
represented in momentum space are employed, namely, 
the Nijm-93, Nijm-I [H, the CD-Bonn ^ and the 
N'^LO [21| potentials, and the Coulomb force is also used 
commonly. In the calculations, the partial waves up to 
Jr < 6 are taken into account. 



In Fig. 21 the hft dependence of calculated ground- 
state energies of ^^O for pi = 12 using the Nijm-93 and 
the CD-Bonn potentials is shown. The unperturbed en- 
ergy which is shown as "unp." and the energy with the 
Iplh correction are displayed separately. The expression 
of the ground-state energy with the Iplh effect is given in 
Eq. IpTzjl as —BE. We see that the effect of the Iplh cor- 
rection has a significant contribution attractively to the 
ground-state energy. If we use the HF wave functions, the 
unperturbed ground-state energies should become more 
attractive. 

We note here that the values of hil at which the en- 
ergy minima are obtained differ from each other between 
the Nijm-93 and the CD-Bonn potentials, and also be- 
tween the unperturbed and the unperturbed plus Iplh 
energies, reflecting differences in the properties of the 
two potentials. In the calculation of ^^O, a value around 
= 14 MeV is often employed as a suitable value of 
hn. This value is very close to that determined by em- 
pirical formula such as hH = A5A~'^/^ -25A~^/^ MeV. In 
the present study, however, we regard the value at which 
the energy minimum is obtained as the optimal one. The 
optimal value should be searched for each state in nuclei. 

Figure illustrates the hfl and pi dependences of the 
groimd-state energy with the Iplh effect for the Nijm-93, 
Nijm-I, the N'^LO, and the CD-Bonn potentials. In prin- 
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FIG. 4: The hQ dependence of calculated ground-state en- 
ergies of for pi — 12 for the Nijm-93 and the CD-Bonn 
potentials. 
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FIG. 5: The ftJl and pi dependences of calculated ground- 
state energies with the Ipl^ eflFect of ^^O for various modern 
nucleon-nucleon interactions. 



ciple, we should take the value of pi as large as possible 
until the results do not depend on pi. When we take as 
pi — 14, the results show fairly good convergence for the 
CD-Bonn potential. As for the Nijm-93, Nijm-I and the 
N'^LO potentials, almost convergent results are obtained 
if we take the value pi = 16. Note that the energy for 
pi = 18 at hn = 14 MeV is calculated for the N^LO po- 
tential in order to confirm the convergence. We see that 
the results for pi = 16 and 18 are almost the same. 

In Table m the final results of the ground-state energy 
with the Iplh effect are tabulated for the four potentials 
together with the experimental value. The binding ener- 
gies per nucleon are also shown. The calculated values 
are for the optimal fiil and pi which can be determined 
from the results as shown in Fig. The results for the 
Nijm 93 and the CD Bonn are the least and most attrac- 
tive, respectively, of the four potentials. This tendency 
can also be observed in the Faddeev-Yakubovsky calcu- 
lations for ^He by Nogga et al . 

It is seen that the calculated ground-state energies are 
less bound than the experimental value. In the present 
calculation, higher-order correlations such as the three- 
body cluster terms have not been evaluated. In addi- 
tion, the real three-body force is not taken into account. 
The inclusion of the real three-body force and the higher- 
order many-body correlations would compensate for the 
discrepancies between the experimental and calculated 
values. Such a study remains as an important task for a 
deeper understanding of nuclear ground-state properties 
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TABLE L The calculated ground-state energies with the Iplh 
effect and the binding energies per nucleon of ^^0. In these 
calculated values, the optimal values of fiil for each interac- 
tion are employed. As for the value of pi, we take as pi = 14 
for the CD Bonn and pi = 16 for the other interactions as 
the sufficiently large value as suggested in Fig. |K1 The exper- 
imental values are taken from Ref. |3^. All energies are in 
MeV. 



16q 


Nijm 93 


Nijm I 


N^LO 


CD Bonn 


Expt. 


-Eg.s. 


-99.69 


-104.25 


-110.00 


-115.61 


-127.62 


BE /A 


6.23 


6.52 


6.88 


7.23 


7.98 



in the present approach. A coupled-cluster calculation 
of the saturation property concerning the binding energy 
and charge radius for ^^O by Mihaila and Heisenberg 
has shown that the calculated result agrees well with the 
experimental value when a genuine three-body force is 
included in the calculation 1321 . 



B. '^N and '^O 

Shown in Fig. is the hfl dependence of calculated 
single-particle energies for the Op states in ^^N and ^^O 
for pi = 12 in the case of the CD-Bonn potential. The 
unperturbed energy and the energy with the lp2h cor- 
rection are displayed separately. The unperturbed single- 
particle energy is given in Eq. H36() and that with the lp2h 
correction is in Eq. (|41|l . We see that the unperturbed 
single-particle energies vary considerably at around the 
typical hil = 14 MeV. However, the single-particle ener- 
gies with the lp2h correction have the saturation points 
at around hfl ~ 14 ~ 15 MeV, depending on the single- 
particle states. Note that the minimum points for the 
ground-state energies of ^^N and ^^O correspond to the 
maximum points of the single-particle energies in Fig. |S1 
It should be remarked that the spin-orbit splittings for 
the Op states in ^^N and ^^O are significantly enlarged by 
taking into account the lp2h correction. This effect has 
already shown in our previous works though the calcula- 
tion was performed perturbatively by taking into account 
second-order diagrams on the isospin basis llEi. 

Figure [7| exhibits the hft and pi dependences of the 
single-particle energies with the lp2h effect. It is seen 
that the pi dependence is weaker than that in Fig. 
This is because the results in Fig.Elare the relative values 
for the binding energies of two nuclei as given in Eq. (|41|l . 
while those in Fig. [S] are not the relative ones. We may 
say that, as for the single-particle energies for the hole 
states, the results for pi = 12 are acceptable as final 
results in the present study. 

In Fig. |S1 the final results of the single-particle ener- 
gies with the lp2h effect for the Op states in '^^N and 
^^O for pi = 12 using the four potentials are shown with 
the values of the spin-orbit splitting energy. The opti- 
mal values of hfl for each interaction are searched for 



TABLE II: The calculated single-particle energies with the 
lp2h effect for pi — 12 in '^^N. The values of the spin-orbit 
splitting energy AEis(Op) = -Es.p. (Opi/2) - Es.p.{0p3/2) are 
also tabulated. In these calculated values, the optimal values 
of hQ for each single-hole state and interaction are employed. 
The results for HQ = 14 MeV are also shown in parentheses. 
The experimental values are taken from Ref. [s^ ■ AH energies 
are in MeV. 

^^N Nijm 93 Nijm I N-'^LO CD Bonn Expt. 
l/2"(0pi/2) -14.21 -14.56 -14.80 -15.71 -12.13 

(-14.12) (-14.56) (-14.80) (-15.73) 
3/2"(0p3/2) -18.71 -19.37 -20.23 -21.63 -18.45 
(-19.26) (-19.86) (-20.17) (-21.22) 
A£;i,(Op) 4.50 4.81 5.43 5.92 6.32 
(5.14) (5.30) (5.37) (5.49) 



the binding energies of "'^^O, "'^^N, and "'^^O in calculat- 
ing the single-particle energies through Eq. (jlTT) . We see 
that the calculated spin-orbit splittings are smaller than 
the experimental values though the differences between 
the calculated and experimental values depend on the 
nucleon-nucleon interactions employed. The magnitudes 
of these discrepancies would be reduced if we include a 
genuine three-body force in the calculation as discussed 
in Kefs. HQ. 

In Tables ^ and IIIII the final results of the single- 
particle energies shown in Fig.|Slare tabulated. The re- 
sults for the typical fifl = 14 MeV are also displayed in 
parentheses for reference. In the case of = 14 MeV, 
we use this value commonly in calculating the binding 
energies of "'^^O, ^^N, and ^^O. It is seen that all the cal- 
culated single-particle energies are more attractive than 
the experimental values. The inclusion of the three-body 
force and the evaluation of higher-order many-body cor- 
relations may compensate for the discrepancies between 
the experimental and calculated values. 

In order to see the accuracy of the calculations, it 
would be worthwhile to apply the present method to the 
few-nucleon systems *He, ■^H and ■^He as similar systems 




ft«(MeV) fto(MeV) 

FIG. 6: The h^l dependence of calculated single-particle en- 
ergies for pi = 12 for hole states in ^®0. The left (right) 
figure is for the proton (neutron) levels which correspond to 
the single-hole states in ^'^N (^''O). The CD-Bonn potential is 
employed. 
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pi dependences of calculated 


single- 



particle energies with the lp2/i effect for the proton (^^N) and 
neutron ('^^O) levels. The CD-Bonn potential is employed. 



TABLE III: Same as Table Ell except for ^^O. 



Nijm 93 Nijm I N'^LO CD Bonn Expt. 



l/2-(0pi/2) 
3/2-(0p3/2) 



-17.52 -17.96 -18.37 -19.34 -15.66 

(-17.51) (-18.00) (-18.34) (-19.27) 

-22.03 -22.80 -23.79 -25.27 -21.84 

(-22.72) (-23.37) (-23.79) (-24.83) 

4.51 4.84 5.42 5.93 6.18 

(5.21) (5.37) (5.45) (5.56) 



to "'^^O, ^^N, and ^^O. As for the few-nucleon systems, 
the binding energies have been calculated precisely by 
various methods |3a |. 

In Fig. El the Kl dependence of calculated ground- 
state energies of "^He, '^H, and "^He is shown for p\ ~ 12 
using the CD-Bonn potential. The expression of the un- 
perturbed ground-state energy is given in Eq. (|35|1 . It is 
noted that the formulae for calculating the unperturbed 
ground-state energies of ''He, '^H and '^He are the same. 
In these cases, only the Osi/2 states of the proton and 
neutron are regarded as the hole states. However, the 
results of the unperturbed energies are different between 
■*He and '^H (^Hc and ''He) because of the A dependence 
of the Hamiltonian as given in Eq. (|30|l . The expres- 
sion of the ground-state energy of ^He with the Iplh 
correction corresponds to Eq. (|37|l . and that for ^H and 
■'He with the lp2h effect is similar to Eq. (|39f) . We see 
that although the calculated results of the unperturbed 
ground-state energies of ^H and ^He are the same, the 
energies with the lp2h effect are different because of the 
charge difference. 

In Table IIVI the calculated ground-state energies of 
^He, ''H, and "'He with the corrections for the optimal 
values of hfl which can be determined from Fig.j^are tab- 
ulated together with the results of the no-core shell model 
(NCSM) [33,|3l| and the experimental values. It has been 
shown that the NCSM results agree well with the results 
obtained by accurate methods for few-nucleon systems 
such as the Faddeev-Yakubovsky calculation [s^. It is 
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FIG. 8: The calculated single-particle energies with the lp2h 
effect for pi = 12 in ^^N and ^^O. The values of the spin- 
orbit splitting are also shown. In these calculated values, the 
optimal values of HQ for each single-hole state and interaction 
are employed. 



seen that our results are less bound by several hundred 
keV than the NCSM results. In the present approach, 
higher-order many-body correlations such as the three- 
body cluster terms are not taken into account. The eval- 
uation of the higher-order many-body correlations would 
gain more binding energy. We may say, however, that 
our result of the charge dependence in the relative en- 
ergy between ^H and "'He is in good agreement with the 
NCSM result and also the experimental value. 

This kind of agreement of charge dependence can also 
be seen in the results for '^^N and ^^O as shown in Ta- 
bles Inland mil The experimental energy difference of the 
ground states between ^^N and '^^O is 3.53 MeV. Our re- 
sult of the energy difference of the single-particle energies 
for the Opi/2 orbits between ^^N and ^^O is 3.63 MeV for 
the CD Bonn. One can see that the results for the other 
potentials also agree well with the experimental value. 
Thus, we may say that the effect of the Coulomb force 
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FIG. 9: The ftfi dependence of calculated ground-state ener- 
gies of "He, ^H, and ^He for pi = 12. The CD-Bonn potential 
is employed. 



FIG. 10: The Ml dependence of calculated single-particle 
energies for pi = 12 in ^'^F and ^''O. The CD-Bonn potential 
is employed. 



TABLE IV: Comparison of the ground-state energies of *He, 
^H, and '^He in the present approximation for pi = 12 with 
those in the NCSM calculations and the experimental values. 
The CD-Bonn potential is commonly used in the calculations. 
The experimental values are taken from Ref. "s^l . All energies 
are in MeV. 
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Expt. 


■'He 
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^H 
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-8.48 


"He 




-26.15 


-26.30 
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£g.s.('He)-i5g 


.B.('H) 


0.75 


0.75 


0.76 



for the pp channel is correctly treated in our particle-basis 
formalism. The Coulomb force effect is also discussed in 
the next subsection for ^^F and ^^O. 
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FIG. 11: The hO, and pi dependences of calculated single- 
particle energies with the 2plh effect in ^^F and ^^O for the 
CD-Bonn potential. 



C. "F and "O 

Figure E| shows the Kl dependence of calculated 
single-particle energies for pi = 12 for the Is and Od 
states in ^^F and ^^O with the CD-Bonn potential. The 
unperturbed energy and the energy with the 2plh cor- 
rection are displayed separately. The definition of the 
single-particle energy with the correction is given in 
Eq. l|in|) . We see that all the unperturbed energies are 
rather unbound and considerably vary at around the typ- 
ical hn = 14 MeV. However, some single-particle states 
become bound at the energy minimum points by taking 
account of the corrections. It should be noted that the 
magnitudes of the spin-orbit splitting with the 2plh ef- 
fect for the Od states are not very different from those 
for the unperturbed part at around Ml — 14 MeV. This 
tendency differs from the case of the deeply bound hole 
states for which the lp2h effect plays an important role 
to enlarge the spin-orbit splittings as shown in Fig. 

In Fig. we show the Ml and pi dependences of 
the single-particle energies with the 2plh effect. We see 
that the pi dependence for the 0^5/2 and lsi/2 states 
shows the good convergence at pi = 12. On the other 
hand, the results for the 0^3/2 states do not necessarily 



converge at pi = 12. Since the 0^3/2 states of the proton 
and neutron are highly unbound, it would be necessary to 
take a larger value of pi in order to obtain the convergent 
results. In the present study, however, we employ the 
values for pi = 12 as the final results of the single-particle 
energies in ^^F and ^^O. 

Shown in Fig. E| are the final results of the single- 
particle energies with the 2plh effect in ^^F and ^^O for 
the four potentials with the values of the spin-orbit split- 
ting energy. The optimal values of Ml are employed for 
the results. On the whole, the calculated spin-orbit en- 
ergies are larger than the experimental values in contrast 
to the hole state case. We may say, however, that the 
calculated results become somewhat smaller if we take a 
larger value of pi, because the Oc?3/2 states are lowered 
as suggested in Fig. ^2 The calculated results for the 
ground Ods /2 states agree fairly well with the experimen- 
tal values. However, in our preliminary estimation, the 
three-body cluster effect for the particle state shows a 
repulsive contribution significantly to the single-particle 
energy, while that for the hole state is essentially small, 
as far as only the two-body interaction is employed. 

In Tables and IVII the final results of the single- 



12 



q- 2 



-2 



17r 



4.84 5.16 5.64 6.35 



3/2+ (Od,,; 



5.00 



■Nijm93 Nijml 

Pi = 12 



CD Bonn 5/2"'(0d5,2) 



17/ 



o 



3/2-'(0Q'3,2r 



5.14 5.57 



6.03 6.78 



5.08 



Nijm 93 



Niim I ^5 



1/2+ (1s„, 



p, = 12 



C^nn 5^2+ (Od,, 



FIG. 12; The calculated single-particle energies with the 2plh 
effect for pi = 12 in ^^F and ^^O for the CD-Bonn potential. 
The values of the spin-orbit splitting are also shown. In these 
calculated values, the optimal values of ftfi for each single- 
particle state and interaction are employed. 



TABLE V: The calculated single-particle energies with the 
2plh effect for pi = 12 in ^^F. The values of the spin-orbit 
splitting energy AEis{Qd) = -Es.p. (Od3/2) - £s.p. (Od5/2) and 
the energy differences between the lsi/2 and Odc,/2 states 
AEsd = -£'s.p.(lsi/2) — ^-s. p. (0^5/2) are also tabulated. In 
these calculated values, the optimal values of hfl for each 
single-particle state and interaction are employed. The re- 
sults for h^l = 14 MeV are also shown in p arentheses. The 
experimental values are taken from Ref. |39|. All energies are 
in MeV. 
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5.14 


5.36 


5.56 
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4.40 


(5.50) 


(5.66) 


(5.43) 
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l/2+(lSi/2) 


0.67 


0.70 


0.58 


0.70 
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(0.84) 
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(0.52) 


(0.54) 
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(0.57) 


(0.80) 




AEisiOd) 


4.84 


5.16 


5.64 


6.35 


5.00 




(5.28) 


(5.45) 


(5.48) 


(5.89) 





TABLE VL Same as TableEl except for ^^O. 



"0 


Nijm 93 


Nijm I 


N^'LO 


CD Bonn 


Expt. 


3/2+(0d3/2) 


2.03 


2.21 


2.30 


2.67 


0.94 


(2.33) 


(2.47) 


(2.19) 


(2.37) 




l/2+(lsi/2) 


-2.55 


-2.57 


-2.82 


-2.76 


-3.27 


(-2.43) 


(-2.49) 


(-2.86) 


(-2.87) 




5/2+(0d5/2) 


-3.11 


-3.36 


-3.73 


-4.11 


-4.14 




(-3.28) 


(-3.33) 


(-3.66) 


(-3.89) 




AE,d 


0.56 


0.79 


0.91 


1.35 


0.87 




(0.85) 


(0.84) 


(0.80) 


(1.02) 




AEisiOd) 


5.14 


5.57 


6.03 


6.78 


5.08 




(5.61) 


(5.80) 


(5.85) 


(6.26) 





particle energies shown in Fig. E| are tabulated together 
with the spin-orbit splitting energies for the Od states 
and the energy differences between the lsi/2 and 0^5/2 
states. The values for the typical hil. = 14 MeV are also 
displayed in parentheses for reference. We may say that 
our results for the magnitudes of these two splittings are 
not very different from the experimental values. How- 
ever, we should take account of the real three-body force 
and evaluate higher-order many-body correction terms 
to obtain more reliable results. This kind of study is in 
progress. 

We here discuss effects of the Coulomb force. The ex- 
perimental mass difference between ^^F and ^^O is 3.54 
MeV. The calculated results lie between 3.41 and 3.73 
MeV for the four potentials in the case of the optimal hfl 
and between 3.50 and 3.63 MeV in the case of hil — 14 
MeV as seen from the values in Tables IVl and IVll These 
calculated values are in good agreement with the exper- 



imental value. This kind of agreement has been shown 
also for the hole states. Furthermore, another effect of 
the Coulomb force appears in the particle states. The 
experimental lsi/2 states of the proton and neutron lie 
above the 0^5/2 states by 0.49 and 0.87 MeV in energy, 
respectively. Thus, the lsi/2 state in ^^F is close to 
the 0^5/2 state by 0.38 MeV than in "O. This effect is 
known as the Thomas-Ehrman shift due to the Coulomb 
force ^1 ■ our results, the magnitudes of 

the shift are from 0.19 to 0.29 MeV for the cases of the op- 
timal hfl and from 0.22 to 0.23 MeV for hQ = 14 MeV, 
depending on the interactions employed. In the latter 
case the results hardly depend on the potentials, because 
the unperturbed lsi/2 h.o. wave functions are the same 
for all the cases using the four interactions, and thus 
the Coulomb force works equally in the calculations. Al- 
though some discrepancies between the experimental and 
calculated values are seen, we may say that the Thomas- 
Ehrman effect can be observed in our results. 
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IV. SUMMARY AND CONCLUSIONS 

The method for calculating the ground-state energy 
and single-particle energy has been developed within 
the framework of the unitary-model-operator approach 
(UMOA). The expressions for the numerical calculation 
have been recast from the isospin basis to the particle one 
for the charge-dependent structure calculation. We have 
applied the UMOA to ^^o, ^^N, ^^O, ^^F, and ^^O em- 
ploying modern nuclcon-nuclcon interactions, such as the 
Nijm-93, Nijm-I, the CD-Bonn, and the N^LO potentials 
which have charge dependence. The Coulomb force has 
been also used for the pp channel. In order to obtain the 
final results, we have searched for the optimal values of 
hfl and the values of pi for which the calculated results 
almost converge. 

The accuracy of the approximation in the present 
method has been investigated by calculating the ground- 
state energies of '^He, '^H, and '^He and comparing the 
present results with the accurate no-core shell-model 
(NCSM) results. We have found that the energy dif- 
ferences between the NCSM and our results for these 
systems are several hundred keV for the CD-Bonn po- 
tential. As for the energy difference between ■^H and 
■^He, our result agrees well with the NCSM result and 
the experimental value. 

The good agreement for charge dependence between 
the present results and the experimental values is ob- 
served also in the differences in the ground-state energies 
between ^^N and ^^O, and and ^'^O. The effect of 
charge dependence is also seen in the Thomas-Ehrman 
shift for the lsi/2 states in ^^F and ^^O. 

We have shown that the calculated spin-orbit splittings 



for the Op hole states are enlarged significantly by taking 
the lp2h effect into account, and become close to the ex- 
perimental value. On the other hand, the influence of the 
inclusion of the 2plh effect on the spin-orbit splittings for 
the Od particle states is rather small. On the whole, the 
calculated spin-orbit splittings for the hole and particle 
states in nuclei around "'^^O are not very different from the 
experimental values though the results somewhat depend 
on the interactions employed. 

In the present work, higher-order many-body correla- 
tions such as the three-body cluster terms are not eval- 
uated. In addition, the real three-body force is not in- 
cluded in the calculations. We should take account of 
these effects for a deeper understanding of the nuclear 
structure. 

By virtue of the extension of the calculation method 
to the particle basis, the present method can be applied 
to proton- or neutron-rich nuclei in the same manner. 
The mechanism of the variation of magic numbers near 
the drip lines may be clarified from a microscopic point 
of view. The study of neutron-rich nuclei around ^^O is 
in progress. Results for these systems will be reported 
elsewhere in the near future. 
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